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ABSTRACT 

We present surface photometry and stellar kinematics of a sample of 5 SBO galax- 
ies: ESO 139-G009, IC 874, NGC 1308, NGC 1440 and NGC 3412. We measured their 
bar pattern speed using the Tremaine- Weinberg method, and derived the ratio, TZ, of 
the corotation radius to the length of the bar semi-major axis. For all the galaxies, TZ 
is consistent with being in the range from 1.0 and 1.4, i.e. that they host fast bars. 
This represents the largest sample of galaxies for which TZ has been measured this way. 
Taking into account the measured distribution of TZ and our measurement uncertain- 
ties, we argue that this is probably the true distribution of TZ. If this is the case, then 
the Tremaine- Weinberg method finds a distribution of TZ which is in agreement with 
that obtained by hydrodynamical simulations. We compared this result with recent 
high-resolution A''-body simulations of bars in cosmologically-motivated dark matter 
halos, and conclude that these bars are not located inside centrally concentrated dark 
matter halos. 
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^ ' 1 INTRODUCTION 



The pattern speed of a bar, Clp, is its main kinematic ob- 
servable. When parameterized by the distance-independent 
ratio TZ = Dl/cib (where Dl is the Lagrangian/corotation 
radius, at which a star is at rest in the bar's rest frame, 
and as is the bar semi-major axis), it permits the classifi- 
cation of bars into fast (1.0 <TZ < 1.4) and slow {TZ > 1.4) 
ones, li TZ < 1.0 orbits are elongated perpendicular to the 
bar, so that self-consistent bars cannot exist in this regime 
(Contopoulos 1980). 

A robust method for measuring TZ relies on hydrody- 
namical simulations to model gas, particularly at shocks. 
These studies find fast bars {e.g. Lindblad et al. 1996; Lind- 
blad & Kristen 1996; Weiner et al. 2001). Hydrodynamical 
simulations can also obtain TZ by matching morphological 
features in H I {e.g. Laine 1996; England et al. 1990; Hunter 
et al. 1989; Aguerri et al. 2001). Moreover, if the leading, 
oiTset dust lanes frequently found in bars can be identified 
with shocks, then fast bars seem to be the norm in late- 
type galaxies (van Albada & Sanders 1982; Athanassoula 
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1992). A model- independent method for measuring flp di- 
rectly was obtained by Tremaine & Weinberg (1984). The 
Tremaine- Weinberg method (hereafter, the TW method) is 
given by the simple expression XQpSini = V, where X and 
V are luminosity-weighted mean position and velocity mea- 
sured along slits parallel to the line-of-nodes. If a number of 
slits at different offsets from the major-axis are obtained for 
a galaxy, then plotting V versus X for the different slits pro- 
duces a straight line with slope flp sin i. To date this method 
has been applied successfully to 3 early-type barred galaxies 
(Kent 1987; Merrifield & Kuijken 1995; Gerssen et al. 1999; 
Debattista et al. 2002) and preliminary results have been 
obtained for another 5 galaxies ranging from SBO to SBb 
(Debattista & Wilhams 2001; Gerssen 2002). Although aU 
these galaxies are consistent with having fast bars, the sam- 
ple size is still small enough that the range of TZ recovered 
by the TW method is still poorly defined. We have started a 
program to enlarge the sample of pattern speeds measured 
with the TW method, including to galaxies of various bar 
strength, environment, luminosity, inclination etc. In a pre- 
vious paper, (Debattista et al. 2002, hereafter Paper I), we 
reported on the special case of NGC 1023, a system which 
shows evidence of a past interaction with one or more of 
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its satellite galaxies. In this paper, we present results for an 
additional 5 SBO galaxies. 



2 THE SAMPLE 

Our preliminary sample of galaxies consisted of 11 SBO ob- 
jects selected to be bright and undisturbed, with no evidence 
of dust or spiral arms to complicate the TW analysis, gen- 
erally not too many nearby stars, and a bar at intermediate 
angle between the orientation of the major and minor axes. 
After detailed analysis of their photometrical properties (as 
described in Section ^ below), 5 of them were rejected due 
to either strong variation in the ellipticity and/or position 
angle of the isophotal profiles at large galactocentric dis- 
tances, or due to the presence of dust in their inner regions. 
Of the remaining 6, Paper I studied one of them, NGC 1023; 
in this paper we will analyze the remaining 5 galaxies: ESO 
139-G009 (SABO), IC 874 (SBO), NGC 1308 (SBO), NGC 
1440 (SBO), and NGC 3412 (SBO). A compilation from the 
literature of their properties is given in Table ^. 

Four of these galaxies, ESO 139-G009, IC 874, NGC 
1308, and NGC 1440 were previously poorly studied. Indeed 
no surface photometry or stellar kinematics was available for 
the first three of these, and only band photometry was 
available for NGC 1440 (Kodaira et al. 1990), from which 
Baggett et al. (1998) obtained a decomposition into bulge 
and disc. 

This was not the case of NGC 3412, for which both 
surface photometry and stellar kinematics were available. 
Optical and NIR imaging of NGC 3412 were obtained by 
Kent (1984, r), Kodaira et al. (1990, V), Shaw et al. (1995, 
JHK) and Ann (2001, UBVRI and Ha). Both Shaw et al. 
(1995) and Ann (2001) measured an isophotal twist in the 
nuclear region of this galaxy, suggesting that its bulge is 
triaxial. Bulge-disc decompositions were perfomed by Kent 
(1985, r) and Baggett et al. (1998, V). Kuijken et al. (1996) 
obtained the position-velocity diagram of the stellar com- 
ponent along PA = 155° (the "major-axis", but see our 
discussion in Section ^ below) of NGC 3412, to search for 
counter-rotation. They found that the fraction of counter- 
rotating disc stars is less than 5 per cent. Spatially resolved 
stellar kinematics of this galaxy have been measured by 
Fisher (1997, PA = 55°, 155°) and Neistein et al. (1999, 
PA = 155°), who derived the TuUy-Fisher relation for a 
sample of nearby lenticular galaxies. 



3 SURFACE PHOTOMETRY 

The photometric observations of ESO 139-G009 and IC 874 
were carried out with the 3.5-m New Technology Telescope 
(NTT) at the European Southern Observatory (ESO) in La 
Silla (Chile) in May 2001. We used the ESO Muhi-Mode In- 
strument (EMMI) in red imaging and low-dispersion spec- 
troscopic (RILD) mode. The detector was the No. 36 Tek- 
tronik TK2048 EB CCD with 2048 x 2048 pixels. The im- 
age scale was 0.27 arcsec pixel"^. We imaged NGC 1308 
and NGC 3412 at the 1-m Jacobus Kapteyn Telescope lo- 
cated at the Roque de los Muchachos Observatory (ORM) 
in La Palma (Spain) in October 2001 and May 2000, respec- 
tively. We used a SITe2 CCD with 2048 x 2048 pixels and 



Table 2. Log of the photometrical observations 



Galaxy 


Date 


Exp. Time (i 


^oc) 






B 


V 


/ 


ESO 139-G009 


22-23 May 2001 


1440" 


480" 


300'' 


IC 874 


22-23 May 2001 


600" 




300'' 


NGC 1308 


11 Oct 2001 






7200= 


NGC 1440 


11 Nov 2001 


2400" 




3600'' 


NGC 3412 


27-28 May 2000 


2400^= 


2700^= 


3600^= 


NOTE - " Bessel, 


' Gunn, Harris 









an image scale of 0.33 arcsec pixel ^ . Finally, we observed 
NGC 1440 with the Danish 1.54-m telescope at ESO in La 
Silla in November 2001. The telescope mounted the Danish 
Faint Object Spectrograph and Camera (DFOSC) with the 
EEV/MAT 44-82 CCD with 2048 x 4096 pixels. The image 
scale was 0.39 arcsec pixel"^. The log of the photometrical 
observations is given in Table ^ 

All images were reduced using standard IRAfQ tasks. 
First, we subtracted a bias frame consisting of 10 exposures 
for each night. The images were then fiat-fielded using sky 
fiats taken at twilight. The sky-background level was re- 
moved by fitting a second order polynomial to the regions 
free of sources in the images. Special care was taken at this 
point to obtain a final background in the images free of 
sources, in order to reach the outermost parts of the ob- 
jects. Cosmic rays were removed by combining the different 
exposures for each filter using a sigma clipping rejection al- 
gorithm. 

Photometric calibrations were achieved by means of 
standard stars. The calibration constants include corrections 
for atmospheric and Galactic extinction, and a colour term. 
Fig. ^ shows calibrated /-band images of the galaxies. 



4 BAR LENGTH 

4.1 Photometrical parameters 

Isophote-fitting with ellipses, after masking foreground stars 
and bad pixels, was carried out using the IRAF task ELLIPSE. 
In aU cases, we first fit ellipses allowing their centers to 
vary; within the errors, no variation in the ellipse centers 
was found for the 5 galaxies studied in this paper. Since 
patchy obscuration can cause variations of ellipse centers, 
we conclude that there is no evidence for such obscuration 
in our sample, as we also verified from colour maps and 
B — I colour profiles, where possible. The final ellipse fits 
were done at fixed ellipse centers. In Figs. we show the 
resulting surface brightness, B — I colour, ellipticity and po- 
sition angle profiles of the galaxies. The inclination and posi- 
tion angle of the galaxies were determined by averaging the 
outer isophotes, as shown in Fig. ^. The values obtained are 
reported in Table ^. The inclinations obtained for the galax- 
ies are, within the errors, consistent with the values reported 
in the literature. The position angle has larger differences, 
the biggest being APA = 99° in the case of NGC 1308. 

1 IRAF is distributed by NOAO, which is operated by AURA Inc., 
under contract with the National Science Foundation 
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Table 1. Parameters of the sample galaxies 



Galaxy Type i PA Bt D25 x das Vcbr D M^^ 

(RC3) (°) {°) (mag) (arcsec) (km s"!) (Mpc) (mag) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


ESQ 139-G009 
IC 874 
NGC 1308 
NGC 1440 
NGC 3412 


(R)SABO(rs) 

SBO(rs) 

SBO(r) 

(R')SBO(rs) 

SBO(s) 


35.6+5-0 
43 6+'^-2 
40 7+''-^ 
55.8_2;8 


94 
17 

135 
28 

155 


14.35 
13.60 
14.70 
12.90 
11.43 


74 X 60 
60 X 43 
70 X 51 
128 X 97 
218 X 123 


5389 
2602 
6180 
1382 
1202 


71.9 
34.7 
82.4 
18.4 
16.0 


-20.28 
-19.41 
-19.88 
-18.90 
-19.73 



NOTE - Col. (2): Morphological classification from RC3. Col. (3): Galaxy inclination from logi?25 in RC3. Col. (4): Major-axis position 
angle from RC3, except for NGC 1308 and NGC 3412 (LEDA). Col.(5): Total observed blue magnitude from LEDA. Col.(6): Apparent 
isophotal diameters measured at a surface- brightness level of = 25 mag arcsec"-^ from RC3. Col. (7): Radial velocity with respect to the 
CMB radiation from RC3, except for NGC 1308 (LEDA). Col.(8): distance obtained as Vcbb./Ho with Hq = 75 km s'^ Mpc"!. Col.(9): 
Absolute total blue magnitude from Mg corrected for inclination and extinction as in LEDA and adopting D. 



ESO 1,"59-G009 IC 874 NGC 1308 





Figure 1. /-band images of the sample galaxies. The isocontours are spaced at 1 mag arcsec"^ and the outermost corresponds to /ij = 23 
mag arcsec"^. In each panel, North is up and East is left, and the bottom horizontal line is 30 arcsec long. For each galaxy, the thick 
continuous line and the dashed lines show the slit position along major-axis and the offsets in the spectroscopic observations. 



We speculate that the LEDA value of the position angle for 
this galaxy is contaminated by the bright star to its east. In 
measuring the photometrical parameters of this galaxy, we 
masked out the whole eastern part, to minimize the contam- 
ination from this star. In the case of NGC 3412, it is possible 
to compare our values of i and PA with the ones obtained 
by Kent (1984), who found i = 56?6 and PA = 151?6. These 
data are in very good agreement with our measurements, as 
given in Tab. |3[ 

Using the ellipse fits, we measured the angle, Vbar, be- 
tween the bar major-axis and the line-of-nodes in the de- 



projected plane of the disc. The values of this parameter are 
reported in Tab. H. 



4.2 Bar semi-major axis length 

Determining the length of a bar is not entirely trivial. This 
is particularly true in early-type galaxies, for which there is 
usually no obvious spiral structure or star formation beyond 
the bar marking its end, so that the luminosity distribution 
of the bar gradually melds into that of the disc. The presence 
of a large bulge may complicate further the measurement of 
as. Fourier decomposition of the galaxy light (Ohta et al 
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Semi — major axis (arcsec) Semi — major axis (oroseo) 



Figure 2. Surface brightness, ellipticity, position angle and B — I colour radial profiles of ESO 139-G009. Continuous, dashed and dotted 
lino refer to 7, V and _B— band data, respectively. The thick lines represent the fits to the /-band ellipticity and position angle of the 
galaxy in the disc region. 



1990; Aguerri et al 2000), identification of a change in the 
slope of the surface-brightness profile along the bar major- 
axis (Chapelon et al. 1992), visual ispection of galaxy im- 
ages (Martin 1995), location of minimum/maximum in the 
ellipticity profile (Wozniak et al. 1995), and modeling the 
structural components contributing to the galaxy surface 
brightness (Prieto et al 2001) are among the several meth- 
ods which have been proposed for measuring as. We have 
used three of these, yielding 4 estimates of as (as.i to a_B,4). 
The mean of these is our best estimate of as and we use the 
largest deviations from the mean for the error estimates. The 



values obtained are given in Table ^. The three methods are 
the following; 

(i) Fourier amplitudes. Since a bar represents a bisym- 
metric departure from axisymmetry, a Fourier decomposi- 
tion is a natural way to measure as ■ The method of Aguerri 
et al. (2000) is based on the ratios of the intensities in the 
bar and the inter-bar region. The azimuthal surface bright- 
ness profile of the deprojected galaxies is decomposed in 
a Fourier series. The bar intensity, /t,, then is defined as: 

= -^0 + -^'2 + ^4 + h (where Jo, I2,Ia and /e are the m = 0, 
2, 4 and 6 terms of the Fourier decomposition, respectively). 
The inter-bar intensity is defined as: Ub — lo — I2 + 14 — le- 
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10 20 30 40 50 
Semi — major axis (arcsec) 



10 20 30 40 50 60 
Semi — major axis (oroseo) 



Figure 3. As in Fig. but for IC 874. 



Table 3. Inclination, disc position angle, and intrinsic bar angle 
relative to the line-of-nodes of the sample galaxies 



Galaxy 



(°) 



n 



'/'bar 

n 



ESQ 139-G009 


46.2 


± 


1.5 


93.5 ± 1.0 


77 


± 


1 


IC 874 


38.7 


± 


1.8 


26.6 ± 2.7 


70 


± 


3 


NGC 1308 


35.9 


± 


2.6 


58.9 ±4.1 


60 


± 


1 


NGC 1440 


37.8 


± 


0.9 


25.9 ± 0.6 


40 


± 


1 


NGC 3412 


55.3 


± 


2.1 


151.0 ±0.9 


84 


± 


2 



The bar region is defined as the region where Ib/hb > 
0.5[max(/b//ib) — min(/6/Jif,)] + min(/i,/Jii,); thus the semi- 



major axis of the bar is identified as the outer radius at 
which Ih/Iib = 0.5[max(/6//ib) — mm{h/Iib)] + rain{Ib/Iib)- 
Athanassoula & Misiriotis (2002) applied this method to 
the analytic models of Athanassoula (1992) (for which as 
are known exactly) and found an accuracy of better than 
8 per cent. In practical applications to early- type galaxies, 
adjustment for the presence of a massive bulge, which in de- 
projection may result in a higher peak in h/hb than that 
due to the bar, may be required {e.g. Paper I). The values 
listed as aB,\ in Table ^ are the ones which were obtained 
using this method. Figure ^ shows relative amplitudes of 
the m = 2,4,6 Fourier components and the bar/interbar 
intensity ratio of the galaxies. 
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Semi — major axis (arcsec) Semi — major axis (oroseo) 



Figure 4. As in Fig. ^ but for NGC 1308. Only /—band data are available for this galaxy and therefore no B — I colour radial profile 
is provided. 



(ii) Fourier and ellipse phases. In their analysis of N- 
body models, Debattista & Sellwood (2000) also made use 
of a Fourier decomposition of the surface density. Here we 
adopt their m = 2 (deprojected) phase method for as, 2, 
namely that a bar can extend only as far out as the phase 
of the m — 2 moment is constant (within the errors, which 
here must include uncertainties in PA and i). In A'^-body 
simulations, secondary structures, such as rings or spirals, 
sometimes lead to an over-estimate of the semi-major axis. 
However, our early-type galaxies were selected to be free of 
spirals and rings. A related, but not equivalent, method uses 
the phases of the ellipse fits, which then gives us estimate 



as, 3 in Table In Figure |7|, we plot the phases of the Fourier 
m — 2 moments and of the ellipses, showing as, 2 and as.s- 

(iii) Decomposition of surface brightness profiles. The sur- 
face brightness profiles of galaxies can be decomposed into 
a number of structural components, each of which may be 
described by some simple analytic model. Barred galaxy 
models are complicated by their non-axisymmetric nature; 
a method for such decompositions has been developed by 
Prieto et al. (2001), which we used here. We assumed four 
different components: bulge, disc, bar and lens. The bulge 
luminosity was modeled with a Sersic's law (Sersic 1968), 
while the disc was modelled by an exponential profile (Free- 
man 1970). Two profile types for the bars were used: el- 
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Semi — major axis (arcsec) 



20 40 60 
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Figure 5. As in Fig. ^ but for NGC 1440. 



liptical (Freeman 1966) and flat (Prieto et al. 1997). The 
best fit for NGC 3412 was obtained with an elliptical bar; 
for the other galaxies, we used flat bars. The lenses were 
characterized by a smooth luminosity gradient with a very 
sharp cut-off (Duval & Athanassoula 1983). The model of 
NGC 1308 did not need to include a lens; for all the other 
galaxies, the fit was improved by including one. 

In all, ten free parameters needed to be fit: the central 
surface brightness and the scale-length of the disc, the effec- 
tive surface brightness, the effective radius, and the Sersic 
index of the bulge, the central surface brightness, the semi- 
major and the semi-minor axes of the bar, the central sur- 
face brightness and the scale- radius of the lens. The flts were 
performed on the major and minor-axes of the bar, after de- 



projecting the image by means of a flux-conserving stretch 
along the disc minor-axis. 

We used an interactive, iterative method, described in Pri- 
eto et al. (2001) to obtain initial estimates of all the param- 
eters. Then, these estimates of the parameters were used 
as input parameters in an automatic fitting routine using a 
Levenberg-Marquardt nonlinear fitting algorithm (Press et 
al. 1992). Figure I shows the best fit of the surface-brightness 
profiles along the major axis of the bar for the galaxies of 
our sample. Table ^ also includes the bar semi-major axis 
measured with this method (03,4). 
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00 



Figure 6. As in Fig. ^ but for NGC 3412. 



Table 4. Bar semi-major axis of the galaxies 



Galaxy aB,i as, 2 a.B,3 a.BA < as > 

(arcsec) (arcsec) (arcsec) (arcsec) (arcsec) 



ESQ 139-G009 


23.4 


14.0 


16.1 


14.4 


-|7 (,+6.4 


IC 874 


18.7 


21.0 


25.0 


15.0 


19.915-1 


NGC 1308 


13.2 


9.0 


13.2 


14.2 




NGC 1440 


23.2 


24.6 


30.5 


19.2 


24.4+^-1 


NGC 3412 


28.2 


30.0 


32.2 


34.0 


31.11^;^ 
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Table 5. Log of the spectroscopic observations 



I -I- a 1 vir 


Date 


Slit 










(sec) 


ESO 139-G009 


24 May 2001 


10 arcsec N 


7200 




Mav 9001 




5400 




23 May 2001 


8 arcsec S 


7200 


IC 874 


23 May 2001 


10 arcsec E 


10800 




22 May 2001 


major axis 


5400 




24 May 2001 


10 arcsec W 


10800 


NGC 1308 


15 Nov 2001 


5 arcsec E 


10800 




17 Nov 2001 


major axis 


5400 




16 Nov 2001 


5 arcsec W 


10800 


iNLrU i44U 


1 / InOV zUOi 


9 arcsec E 


7200 




15 Nov 2001 


6 arcsec E 


5400 




15 Nov 2001 


major axis 


5400 




16 Nov 2001 


6 arcsec W 


5400 




17 Nov 2001 


11 arcsec W 


7200 


NGC 3412 


01 Feb 2001 


10 arcsec E 


7200 




02 Feb 2001 


major axis 


7200 




01 Feb 2001 


10 arcsec W 


7200 




02 Feb 2001 


15 arcsec W 


7200 



5 LONG-SLIT SPECTROSCOPY 

The spectroscopic observations of ESO 139-G009, IC 874, 
NGC 1308 and NGC 1440 were carried out with NTT 
at ESO in La Silla in May and November 2001. The 
NTT mounted EMMI in red medium-dispersion spectro- 
scopic (REMD) mode, using the grating No. 6 with 1200 
grooves mm~^ in first order with a 1.0 arcsec x 5.5 arcmin 
slit. The detector was the No. 36 Tektronix TK2048 EB 
CCD with 2048 x 2048 pixels of 24 x 24 ^m^. It yielded 
a wavelength coverage between about 4840 A and 5490 A 
with a reciprocal dispersion of 0.320 Apixel"^. The instru- 
mental resolution was 1.19 A (FWHM) corresponding to 
o'inst ~ 30 km at 5170 A. The spatial scale was 0.270 
arcsec pixel" ^. 

We observed NGC 3412 with the 3.6-m Telescopio 
Nazionale Galileo (TNG) at the ORM in La Palma in Febru- 
ary 2001. The TNG was equipped with the Low Resolution 
Spectrograph (DOLORES); we used the HR-V grism No. 6 
with 600 grooves mm~^ in combination with the 0.7 arcsec 
X 8.1 arcmin slit and the Loral CCD with 2048 x 2048 pix- 
els of 15 X 15 fim^. The wavelength range between about 
4700 and 6840 A was covered with a reciprocal dispersion of 
1.055 A pixeP^. The instrumental resolution (obtained by 
measuring the width of emission lines of a comparison spec- 
trum after the wavelength calibration) was 3.10 A (FWHM). 
This corresponds to an instrumental dispersion of ai„st ~ 80 
km at 5170 A. The spatial scale along the slit was 0.275 
arcsec pixel" ^. 

For each galaxy, we took a spectrum with the slit along 
the major axis, as well as two or more offset spectra with 
the slit parallel to the major axis on each side of the center. 
The major axis position angles were chosen according to the 
values obtained from the surface photometry, as described 
in Sect, k.ll and listed in Table H. The location of the off- 



set spectra was chosen, using the photometry, to give large 
X (and therefore large V) and high S/N. Integration time 
of the galax;y spectra was split into exposures of 2700/3600 
sec to deal with cosmic rays. A log of the observations is 
given in Table ^. In each observing run we obtained spectra 
of twelve giant stars with spectral type ranging from late- 
G to early-K, selected from Faber et al. (1985), to use as 
templates in measuring stellar kinematics. Arc lamp spec- 
tra were taken before and/or after every object exposure to 
allow an accurate wavelength calibration. The range of the 
seeing FWHM during the observing runs was 0.7-1.6 arcsec 
in February 2001 (as measured by fitting a two-dimensional 
Gaussian to the TNG guide star), 0.5-2.2 arcsec in May 
2001, and 0.5-1.5 arcsec in November 2001 (as measured by 
the ESO Differential Image Motion Monitor). 

All the spectra were bias subtracted, Hatfield corrected, 
cleaned of cosmic rays, corrected for bad pixels and columns, 
and wavelength calibrated using standard MIDAS^ routines as 
in Paper I. We checked that the wavelength rebinning was 
done properly by measuring the difference between the mea- 
sured and predicted wavelengths (Osterbrock et al. 1996) 
for the brightest night-sky emission lines in the observed 
spectral ranges. The resulting accuracy in the wavelength 
calibration is better than 2 km . 

The spectra taken along the same axis for the same 
galaxy were co-added using the center of the stellar contin- 
uum as reference. The contribution of the sky was deter- 
mined from the outermost ~ 30 arcsec at the two edges of 
the resulting spectra, where the galaxy light was negligible, 
and then subtracted, giving a sky subtraction better than 
1 per cent. A one- dimensional sky-subtracted spectrum was 
obtained for each kinematical template star. 



6 BAR AND DISC KINEMATICS 
6.1 Stellar kinematics 

We measured the stellar kinematics from the galaxy absorp- 
tion features present in the wavelength range centered on 
the Mg line triplet (AA 5164, 5173, 5184 A, see Tab. us- 
ing the Fourier Correlation Quotient method (Bender 1990; 
Bender et al. 1994), as done in Paper I. We adopted HR 6817 
(KlIII) as the kinematical template to measure the stellar 
kinematics of NGC 3412, HR 7429 (K3III) for the kinemat- 
ics of ESO 139-G009 and IC 874, and HR 3145 (K2III) for 
the kinematics of NGC 1308 and NGC 1440. 

The values of line-of-sight velocity v, and velocity dis- 
persion a measured along the different slits for each sample 
galaxy are given in Table H. 



^ MIDAS is developed and maintained by the European Southern 
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Figure 7. Bar semi-major axis length of the sample galaxies. Left panels (from top to bottom): Relative amplitudes of the m = 2,4, 6 
Fourier components, bar/intcrbar intensity ratio, phase angle of the m = 2 Fourier component, and position angle of the dcprojectcd 
isophotal ellipses. The dashed vertical lines show the bar length (as.i, ag i, and 3 of Table 4). Right panels (from top to bottom): 
decomposition of the radial surface brightness profiles in /-band along the bar major-axis, and residuals of the fit. The adopted structural 
components are bulge (dashed line), disc (dotted line), bar (dash-dotted line) and lens (long-dashed line). The continuous line represents 
the total model. 
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Figure 8. The major-axis radial profiles of the stellar line-of-sight velocity dispersion and velocity (after subtraction of the systemic 
velocity) of the sample galaxies. The profiles are folded around the center. Filled dots and asterisks refer to receding and approaching side 
of the galaxy, respectively. For each galaxy the exponential radial profile fitting the line-of-sight velocity dispersion profile and circular 
velocity obtained by applying asymmetric drift correction are plotted in upper and lower panels, respectively. Note that all data at 
\cp\ < 30° from all slits have been used for the fits to V^,flati but only points on the major-axis are shown in this figure. 
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Table 6. Stellar kinematics of the sample galaxies. 
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6.2 Pattern speed measurement 



6.4 Asymmetric drift correction 



We measured the bar pattern speed for each sample galaxy 
by applying the TW method as done in Paper I. 

To compute the mean position of stars, X, along the 
slits, we generally extracted profiles from / and V-band 
surface photometry along the positions of the slits. In the 
case of NGC 1440, because of the radial colour gradient, we 
preferred to use the intensity from the slit spectra them- 
selves, although they are somewhat noisier. For the remain- 
ing galaxies, the /-band profiles match very well the profiles 
obtained by collapsing the spectra along the wavelength di- 
rection, confirming that the slits were placed as intended. 
We used the broad band profiles to compute X because 
these are less noisy than the spectral profiles, particularly 
at large radii. We computed the value of X at each slit posi- 
tion by Monte Carlo simulation, with photon, read-out and 
sky noise to compute the errors. Formally, the integrals in 
TW equation are over — oo < X < oo, but can be limited to 
—Xmax < X < Xmax if Xmax has reached the axisymmetric 
part of the disc; still larger values of Xmax add noise only. 
The X values thus obtained are given in Tab. ^ 

To measure the luminosity-weighted line-of-sight stel- 
lar velocity, V, for each slit position, we collapsed each two- 
dimensional spectrum along its spatial direction to obtain a 
one-dimensional spectrum. The resulting spectra have been 
analysed with the FCQ method using the same template 
star adopted in Sect. pTl to derive the stellar kinematics 
of the galaxy; V is then the radial velocity derived from 
the LOSVD of the one- dimensional spectra. For each slit 
position, the uncertainty on V was estimated by means of 
Monte Carlo simulations as done for v in measuring the stel- 
lar kinematics. The V values we derived along each slit and 
the adopted wavelength (always including the Mg triplet) 
and radial ranges (limited by the noise, and after removing 
the contribution of foreground stars by linear interpolation) 
are given in Tab. ^. 

For each galaxy, we derived fip sin i by fitting a straight 
line to the values {X,V) obtained from the available slit po- 
sitions (Fig. ^). Finally the value of Q,p (Table ^) was derived 
from the adopted galaxy inclination, as given in Table ^. 



Measurement of Q,p with the TW method requires no mod- 
eling. However, in the absence of gas velocities at large radii, 
measurement of TZ requires some modeling to recover the ro- 
tation curve from the observed stellar streaming velocities. 
This asymmetric drift correction can be fairly large in early- 
type disc galaxies, where the velocity dispersions are large. 
We start from the asymmetric drift equation {e.g. Binney & 
Tremaine [1987] Eqn. 4-33) 



2 



9 In p d In (7^ 
dlnR ^ aini? 



(1) 



where Vc is the circular velocity, p is the disk's volume den- 
sity, and and ajj are the tangential and radial velocity 
dispersions in the cylindrical coordinates of the galaxy's in- 
trinsic plane. We then make the following assumptions: 

(i) We assume that the true rotation curve is flat at 
large radii. This is usually a very good approximation 
in high surface-brightness galaxies. Thus we only measure 
Vc,flat, the amplitude of the rotation curve. A flat rotation 
curve implies o-'^/i^r = |, which allows us to write, for 
the observed velocity dispersion on the major-axis, (Tobs = 
o",/, sin iy/l + Ic? cot^ i, where a — az/cFR. We use 3 values 
of a, 0.7, 0.85 and 1.0, which bracket a plausible range for 
this parameter in early-type disc galaxies extrapolating from 
the only 3 measured values in later-type systems (Dehnen 
& Binney 1998; Gerssen et al. 1997; Gerssen et al. 2000). 
The variations in this parameter are included in our error 
estimate of 14, flat- 

(ii) We assume that aobs decreases exponentially with ra- 
dius, with scale-length _Ro-. We flt with the major-axis 
data, outside the bar region when the S/N permits; these 
flts are shown in Fig. |^. 

(iii) We assume the the disc's volume density is radially 
exponential. Since we are interested here in the total mass 
density, we fltted an exponential to the surface brightness 
measured in ellipses with position angle and ellipticity fixed 
to the values determined from large radii. The fits to 
started outside as; we varied the radial range of the fits to 
obtain error estimates on Rd\ these were as large as 15 per 
cent. Variations in i?d are also included in our error estimate 

on Vc,flat- 



6.3 Systemic and mean streaming velocities 

To measure the systemic velocity of the galaxies, we fit inter- 
polating splines to the measured velocities and fitted "tilted 
rings", of fixed i and PA (as obtained from the photome- 
try), to the inner parts of the galaxies. In the case of NGC 
1308, our data were not of sufficient S/N to permit such fits, 
in which case we used V on the major axis; for the other 
four galaxies, comparision of Tables and |^ shows that the 
major- axis V is a very good approximation to Vsya. In Table 
^ we compare our measurements of Vsys with those in RC3. 
In all cases we find excellent agreement. 

Once Vsys is measured, we obtain the stellar mean 
streaming velocities, 14, by subtracting Kys, then folding 
about the origin the major-axis data. Figure ^ shows that 
our folded spectra have no substantial asymmetries on the 
two sides. 



With these assumptions, the asymmetric drift equation 
becomes 



J2 

"obs 



sin^ z(l + cot^ i) 



2r(^— + 



1 



2 

Ra 



(2) 



We apply this correction to all velocity data points (in- 
cluding those not on the major-axis) outside the bar radius 
and within 30° of the major-axis. We average over all these 
points to obtain 14, flat, with an additional error from the 
scatter of the points added in quadrature. We report our 
findings of Vc.flat in Table ^, and in Fig. || 

In Paper I, the high surface brightness of NGC 1023 
allowed us to drop the assumption that (t^/o-r = i, and 
instead we used the velocity dispersion data from all slits to 
obtain 14, flat ~ 270±31 km (having assumed 0.0 < a < 
1.0; because of its high inclination 14, flat for this galaxy is not 
very sensitive to the value of a). For comparison, we have 
redone the asymmetric drift correction of NGC 1023 with 
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Figure 9. Pattern speed measurements for our sample of galaxies. The straight line is the best fit, with slope f2psini. The resulting 
values of flp are given in Table H 
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Table 7. The values of V and X for each sUt of our sample of galaxies. 



Galaxy 


Offset 


X 


V 


1 1 max 


-^min ! -^max 






(arcsec) 


(km s~^) 


(arcsec) 


(A) 




10 arcsec In 




04UD.4 ± ( .u 


A O 
4U 


OU14. i,04oD.z 




major axis 


-0.12 ±0.02 


5375.0 + 6.2 


60 


5014.1,5486.2 




8 arcsec S 


-1.62 ±0.08 


5366.0 + 8.4 


50 


5014.1,5486.2 




10 arcsec E 


— l.oD it U.U4 


ZoUO.D it 4.0 


OO 


OUoy.Z,04oD.Z 




major axis 


ifi 01-1-0 fif^ 




SO 


OUOy .Z,0^c5D.Z 




10 arcsec W 


1 o n7 -U o OR 
-\-z.yj ( ± u.UD 


O'iO'i 1 _l_ /I o 
ZoZo. i It 4.U 




OUoy .Z,04oD.Z 


NGC 1308 


5 arcsec E 


-1.51 ±0.01 


6358.1 + 14.7 


35 


5115.3,5431.7 




major axis 


— U.Ui it U.UU 


bozo. ( it lo. ( 


OO 


Oli0.r5,04oi. / 




5 arcsec 


H-i.40 it U.UI 


DZoo.D it io. 1 


OO 


OiiO.O,04v5i. / 


NOG 1440 


9 arcsec E 


-2.35 ±0.10 


1625.1 + 5.7 


50 


5115.3,5431.7 




6 arcsec E 


-1.46 ±0.10 


1603.5 + 3.5 


50 


5115.3,5431.7 




major a^xis 


+0.25 ±0.10 


1599.8 + 3.2 


50 


5115.3,5431.7 




6 arcsec W 


+2.21 ±0.10 


1594.8 + 3.9 


50 


5115.3,5431.7 




11 arcsec W 


+4.47 + 0.10 


1581.2 + 7.9 


50 


5115.3,5431.7 


NGG 3412 


10 arcsec E 


-1.46 + 0.13 


843.2 + 1.8 


80 


4964.2,5541.4 




major axis 


-0.06 ±0.03 


849.1 + 1.0 


80 


4964.2,5541.4 




10 arcsec W 


+1.19 + 0.07 


852.3 + 2.0 


80 


4964.2,5541.4 




15 arcsec W 


+0.83 + 0.06 


853.5 + 4.0 


80 


4964.2,5541.4 



the method described above for the new sample. We found 
Vc,flat = 270 ± 13 km , which is in excellent agreement 
with the value of Paper I. 

We finally measured TZ as 14,flat/(asf^p)- We used 
Monte Carlo to estimate the uncertainties in TZ by vary- 
ing 14, flat and as uniformly in their respective ranges, and 
varying ftp assuming the errors of Table ^ are Gaussian. We 
report the median and the 67 per cent interval of TZ in Table 



7 DISCUSSION AND CONCLUSIONS 

The 5 SBO's presented in this work, together with NGC 
1023 studied in Paper I, which we include in our sample 
for this discussion, represent the largest sample of barred 
galaxies, with Qp measured by means of the TW method. 
For all of them, TZ is consistent with being in the range 
1.0 to 1.4, within the errors, i.e. with each having a fast 
bar. The unweighted average for the sample is 7?. = 1.1. The 
apparent range of TZ spans from 0.8 to 1.6 (0.6 to 2.1, within 
the 67 per cent intervals). This spread is not related to the 
properties of the galaxies in any obvious way {e.g. for the 
two galaxies at Vc,flat = 277 km NGC 1023 and NGC 
1440, the measured values of TZ are at opposite extremes of 
the distribution). The fact that some of the values of TZ are 
nominally less than unity leads us to suggest that the large 
range of 7?. is a result of random errors and/or scatter in the 
measurements. 

The sources of random errors are largely due to mea- 
surement uncertainties in all 3 quantities used to compute 
TZ, i.e. Vc,flat, as and Qp. Of these, the largest is in Qp, 
amounting to typical fractional uncertainties of 30 per cent, 
followed by as, for which the typical fractional uncertainty 
is 20 per cent. These uncertainties account for the typical 
large (and asymmetric towards large values) errors on mea- 
surements of TZ. 




Figure 10. The corotation radius, Dj^, and the bar semi-major 
axis, as, for the sample galaxies, including NGC 1023 (Paper 
I). The open circles represent galaxies which are not part of our 
sample: NGG 936 (Merrifield & Kuijken 1995) and NGG 4596 
(Gerssen et al. 1999). Dashed lines corresponding to 7?. = 1 and 
TZ = 1.4, separate the fast-bar, slow-bar and forbidden regimes. 
The dotted line corresponds to 7?. = 1.7. 



A likely source of scatter is errors in the disc PA. De- 
battista (2002) shows that, for PA errors of FWHM 5° (note 
that the root-mean-square PA uncertainty in our sample is 
2?1), the scatter in TZ is of order 0.44, large enough to ac- 
count for the TZ < 1 cases. Since PA errors scatter TZ to both 
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Table 8. Disc and bar kinematics 



Galaxy Vsys V^ys^Rcs K.flat -Dl 

(hm s~^) (km s~^) (km s^^) (km arcsec~^) (arcsec) 

ESQ 139-G009 5377 ± 6 5370 ± 60 3141^2 21.4 ±5.8 14.7t3^ ^-^to.l 

IC 874 2309 ±4 2311 ± 36 187+^3 7.0 ± 2.4 26.9+g" ^-^tol 

NGC 1308 6329 ± 16 6350 ± 42 347^18 39.7 ±13.9 8.9l2:3 ^-^-o.t 

NGC 1440 1601 ±2 1534 ± 66 283^" 7.4 ± 1.7 38.31^°/ ^-^-o.l 

NGC 3412 850 ± 2 865 ± 27 205^1® 4.4 ± 1.2 46.91^9^ ^-^-0.3 



larger and smaller values, then the largest measured value 
ofTZ = 1.53 is probably an over-estimate. If this is the case, 
then TW measurements are finding the same range of TZ as 
do hydrodynamic simulations. 

The conclusion that bars are fast constrains the dark 
matter distribution in disk galaxies. Debattista &i Sellwood 
(1998, 2000) argued that bars this fast can only survive if 
the disc in which they formed is maximal. Recent high res- 
olution A'^-body simulations with cosmologically-motivated 
dark matter halos produce bars with TZ in the range between 
1.2 and 1.7 (Valenzuela & Klypin 2002). Even discounting 
our argument above in favor of a more restricted range of 
TZ, Fig. ^ shows that TZ = 1.7 is possible only for the bars 
of IC 874, NGC 1440 NGC 3412 and, marginally, NGC 936, 
while the bars of ESQ 139-G009, NGC 1023, NGC 1308 and 
NGC 4596 never reach this value of TZ. Note, moreover, that 
3 of the galaxies that do reach 7?. = 1.7 have amongst the 
largest fractional errors in TZ. Therefore we conclude that 
the A'^-body models of Valenzuela & Klypin (2002) probably 
produce slower bars than the observed. 

The galaxy ESQ 139-G009 is classified as SAB {i.e. 
weakly barred) in RC3; the fact that it hosts a fast bar 
suggests that weak bars form via the same mechanism that 
forms the strong ones. Thus the hypothesis of Kormendy 
(1979), that weak bars are the end state of slowed down fast 
bars, already in question from A'^-body simulations (Debat- 
tista & Sellwood 2000), is also unsupported by the limited 
observational data. Further measurements of TZ for weak 
bars would be of considerable interest. 
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